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A USE OF LIGHTHILL'S GENERALIZED FUNCTIONS IN THE 
SOLUTION OF OPTIMAL CONTROL PROBLEM 

Ernest S. Armstrong 

The usual definition of the Dirac "6- function" as a function which 
vanishes everywhere except at a single point but whose integral over the 
entire real line is unity is not a rigorous mathematical definition. No 
mathematical function can have these properties. 

There are, however, several methods of rigorously using such symbolic 
expressions. In each the concept of a function is generalized to include 
all continuous functions, all Lebesque integrable functions, and new 
objects a simple example of which is the "6- function" . As a result of 
this generality the theory is referred to as the theory of generalized 
functions. 

Probably the most mathematically sophisticated approach, requiring a 
knowledge of measure theory and Lebesque integration, is due to L. Schwartz. 
An introduction to Schwartz fe French work (reference l) has been prepared in 
English by I. Halperin (reference 2). Fortunately, there is another for- 
mulation, presented in a text by M. T. Lighthill (reference 3) and 
originally due to G. Temple (reference 4), requiring only a knowledge of 
ordinary calculus. Lighthill' s approach, a brief introduction to which 
is presented in Appendix A, is to construct generalized functions in 
terms of sequences of continuous functions. For example, a representation 
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(not unique) for the Dirac "6- function" is the sequence {e- nt2 (V")*}. 

oo 2 1 +2 i 

In the limit as n -» oo ^ e“ n ^ (n/n) 2 dt -+ 1 and e -n (n/n)^ takes on 

-OO 

the properties of the function" . In Lighthill's text it is 

essentially established that sets of equivalent sequences (see Appendix A) 

form an additive abelian group over which differentiation is defined. It 

was found that this theory could be applied to a class of optimal control 

problems. Appendix A contains the basic definitions and properties of 

Lighthill's formulation which are needed to present the application. 

The general result of applying an optimal control theory, such as 

Pontryagin's maximum principle (reference 5), to a controllable dynamic 

system is a two-point boundary value problem. A large class of these 

boundary value problems take the form. 

I. Given the dynamic equation (written in scalar form for 
simplicity) 


x = f(x, sgn p, a, t), t€ [t Q , tj , x(t Q ) = x°, 
determine the parameter at such that at t^, e(t]_) = 0. 


The variable x is a state or conjugate variable of the control theory, 
p is a continuous switching function, and p is a function of x(t^), 
a, and t^ representing a terminal condition. The function f is 
continuous in x, sgn p (see Appendix A for definition), and u, and 
piecewise continuous in t with points of discontinuity occurring at the 
zeros of the continuous function pix(t), u, t'i , The usual approach for 
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the solution of such problems is to devise an algorithm whereby an 

assumed value of a can be corrected so as to drive e(t^) to zero. 

This requires the evaluation of the expression — (t-i) describing the 

da 

influence of a on x(t^). It is the purpose of this discussion to 

indicate how Ligh thill's generalized function theory can be used to 

obtain an equation for — (t), t€ t Q , t-, 1 . 

da L 

In the sense of ordinary functions we can write 


4£ (x, 
da 

Let t* 
p on 


sgn p, a, t) = 

be a zero of 
[to. tj . If 


Of d(sgn p) dp 4. of dx + 9f 

6 (sgn p) dp da Sx da Qa 

p(x(t), a, t) and S(t*) be the set of zeros of 

£P Vf 0 then as a generalized function 
dt t=t* 


££ d(sgn p) dp j_s represented 

5 (sgn p) dp da 


by I Q(t*) t>( t - t*) where 
S(t*) 

dp 

Q(t*) = - 2 sgn p (t*“) Of (t*) ^ (t*) 

5(sgn P } 

6 ( t — t*) = the Dirac "6- function" of t - t* 

and I represents the sum over all t* of p(x(t), a, t) in It , ti 
d( t») L 
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Let P(t) be the generalized function corresponding to the ordinary- 

function 6f dx + Of . As generalized functions we can write 

6x da 5a 

by property 3 of Appendix A. 


$£=L(x)=l Q(t*) &(t - t*) + P(t). 
dt da da S(t*) 


Applying first property 1 and then 2 of Appendix A the (ordinary function) 
% 

equation results . 


5* 1 

t„ 


r? ^ (t 0 ) - 1 Q(t*) H(t - 1^) 
L da m s(t*) 


t 0f 

5 (1 


i * g.] 


dF 

dt 


dt = 0 


for all "good" functions F(t). 


(see Appendix) 

This in turn leads, by the integral form of the fundament lemma of the 
calculus of variations (reference 6) , to the integral equation 


II. 


£ (t) = dx (t ) 
da da 


+ 1 Q(t») H(t - t*) + 
S(t») 


r /of dx . 5f . 

{ si )d3 - 


By solving this equation in conjunction 


dx (t]_) results, 
da 


with the x equation, 
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By way of application of generalized functions in a computational 

algorithm, the following problem was solved: 

find the optimal thrust ma gn itude and direction of 
a space vehicle, launched from the surface of the 
moon, such that it will rendezvous with a minimum 
of fuel expenditure with a target vehicle in a 
circular orbit about the moon. The moon is assumed 
to have an inverse square law gravitational field. 

Rendezvous is interpreted to be the nulling of the 
relative position and velocity of the two vehicles. 

Pontryagin's maximum principle (reference 5 ) was applied to the 

three dimensional equations describing the relative motion of the two 

spacecraft yielding an optimal control law and a (vector) two point 

boundary problem of the form I. An algorithm was formulated requiring 

the evaluation of dx(t^) (now a matrix) which was computed by an 

da 

integral equation of the form II. The equations were programed and 
solved on an IBM 7094 digital computer. The algorithm was shown to work 
quite well yielding both planar and non -planar minimum fuel rendezvous 
trajectories in about 7 minutes of computer time. 

Space does not permit a complete discussion of the algorithm and 
optimization results. It is intended that this discussion will serve 
to make the reader aware of a particular use of Ligh thill's theory of 
generalized functions in analysis and of a specific application resulting 
from the use of optimal control theory. More complete details of both the 
theory and application are found in reference 7. 
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APPENDIX A 


GENERALIZED FUNCTIONS 


Certain definitions and properties of generalized functions as 
presented in reference 3 are contained in this section for reference 
in the main text and to give the reader a general idea of the con- 
cepts involved. 


Definition 1 - 


A good function. F(t), is one which is 
everywhere differentiable any number of 
times such that it and all its derivatives 
approach zero as its argument approaches 

+oo. 


Definition 2 - A sequence h n (t) of good functions is 

called regular if, for any good function 
F(t), the limit 


lim 


J h n (t) F(t) exists. 


Definition 3 ~ 


A generalized function, h(t), 
as a regular sequence h n (t) of go 
functions, i.e. h(t)-» [h n (t)J • 


is defined 
good 


Thus each generalized function is really 
the class of all regular sequences equivalent, 
in the sense that the limit in Definition 2 
is the same for each sequence, to a given 
regular sequence. 


By the symbol j°° h(t) F(t) dt one means 


^ f hn(t) F (t) dt. 
n-*°° -so 
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Definition 4 - Two generalized functions, h(t) and 

g(t), are equal if and only if 

f° h(t) F(t) dt = j°° g(t) F(t) dt 

-CO Zco 

for all good functions F(t). 

Definition 5 - The derivative h* (t) is defined by the 

sequence h' n (t), i.e. h* (t)-» [V (t)] 


Property 1 h’ (t) F(t) dt = - f° h(t) F' (t) dt 

—CO 

Property 2 If the ordinary function f (t) is such that 

j°° (t^+l)” n f(t) dt exists (for sane natural 
number N) then a generalized function h(t) 
exists such that for all F(t) 


r h(t) F(t) dt = r f(t) F(t) dt. 

-CO -CO 


If sgn(t) 


The integral on the right is the integral in the ordinary 
sense. When the generalized function h(t) has been 
defined, this integral has a meaning also in the theory 
of generalized functions and the above equation states that 
these two meanings are the same. 


r l , t >o 
l-i , t<o 


Signum function 


V 


H(t) = 1/2 (l + sgn t) 


Heaviside function 
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0 , t#0 

at t=0, 6(t) is infinite in such 
/ a way that 

i r 6 (t) dt = 1 

V -oo 


Dirac "6- 


then d sgn t = 2 dfl = 26 (t) 
dt dt 

Definition 6 - If h a (t) is a generalized function of 

for each value of a parameter a then 


cMi^t) = lim ha 1 (t) - hq(t) 
da a'-*a a' - a 


Property 3 If dhq(t) exists, d_ d hq ( t ) 

da da dt 

= L. ±_ h a (t) 
dt da 


function" 



